Optimization of spin-torque switching using AC and DC pulses 
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We explore spin-torque induced magnetic reversal in magnetic tunnel junctions using combined 
AC and DC spin-current pulses. Employing a macrospin model we calculate the optimal pulse times 
for both the AC and DC pulses, needed to minimize the Joule heat loss during the switching process 
for a given AC driving frequency and spin-current magnitude. The results of this optimization are 
compared against numeric simulations. Finally we show this optimization leads to different dynamic 
regimes, where switching is optimized by either a purely AC or DC spin-current, or a combination 
AC/DC spin-current, depending on the anisotropy energies and the spin-current polarization. 



PACS numbers: 85.75.Dd, 85.75.-d, 75.75.Jn 



The question of spin-torque (ST) efficiency in switch- 
ing the free layer in magnetic tunnel junctions (MTJ) has 
been of great interest in recent years [ll41| due to the po- 
tential of spintronic memory as a universal non- volatile 
memory element [1, [H]- For co-linear MTJs, using DC 
current pulses, the most efficient prescription has been 
shown to be roughly twice the critical current [1, 0) ■ Re- 
cent experimental research, however, has shown perpen- 
dicularly polarized DC spin-current pulses can produce 
switching much faster and more efficient than in the co- 
linear case [1,0. Furthermore, previous theoretical work 
by the authors of this paper has shown perpendicularly 
polarized AC spin-currents can improve the efficiency fur- 
ther in MTJs with strong easy-plane anisotropy by induc- 
ing a resonant response in the free layer (lol |. 

For MTJs with weak easy-plane anisotropy (such as 
those used in Ref. pll - fisj ) however, this purely AC 
method becomes less effective. This is a result of the 
magnetization spending more time at large azimuthal 
angle through the switching process where the perpen- 
dicular ST is weaker. Conversely, as the magnetization 
spends more time at large azimuthal angles the strength 
of the parallel DC ST increases. In fact, for some cases 
the strengths of the AC and DC ST may intersect allow- 
ing AC ST to dominate for low energy orbits and DC ST 
for high energy orbits. This tradeoff suggests an alter- 
native means of magnetic switching in MTJs with weak 
easy-plane anisotropy. Instead of using a purely AC or 
DC spin-current, an AC pulse may be used to push the 
magnetization to a higher energy state, where a DC spin- 
current can then be used to switch the magnetization the 
rest of the way. Such an AC/DC current pulse strategy is 
considered experimentally in Ref. and theoretically 
using micromagnetic simulations in Ref. [l5| and shown 
to markedly improve the efficiency of the switching pro- 
cess. 

In this paper provide a theoretical description of mag- 
netic switching using such AC/DC spin-currents. This 
description is then use to determine the optimal AC 
and DC current pulse times which minimize the Joule 
heat loss during the switching process. Such switching is 
then optimized numerically with respect to spin-current 



strength and AC drive frequency which is then compared 
to numeric simulations. The macrospin model employed 
is a simplification of the free layer dynamics. It is useful 
to provide intuition and approximate a semi-quantitative 
description of the switching kinetics. Good agreement 
between macrospin models and experimental results for 
AC and DC switching strategies are found in Ref. [T^ . 
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FIG. 1: (Color online) Simulated optimal Joule heat loss 
during switching process in arbitrary units as a function 
of spin current amplitude (Xs/Ms) for DC (black squares) 
and AC/DC (red triangles) spin current methods. Black 
line represents the theoretical optimal Joule heat loss. Here 
a = 0.015, = 0.5Ms, = 0, if'^ = Jf"^ = Is, 
nip — —Bx — and T = SOOiS'. Inset shows diagram of 
MTJ. Blue and red regions represent fixed and free layers 
respectively. Yellow regions non-magnetic and current flows 
along vertical axis. Arrows indicate fixed layer polarizations 
and free layer easy axis. 

To model the magnetic switching using this AC/DC 
strategy we treat the free layer as a single magnetic do- 
main with a constant saturation magnetization and 
magnetization direction specified by a time-dependent 
unit vector m(<). Its dynamics is described by the 
Landau-Lifshitz-Gilbert equation [l6l | with Slonczewski 
spin torque term p7l . [l8| . Thermal noise is included 



2 



as a random contribution h(t) to the effective field 
Hefi(m, oc —ViiiE{m) with correlator given by 
the the fluctuation-dissipation theorem [23|- Here E is 
the magnetic anisotropy energy 



,(1) 



where /xq is the permeability of free space, and e^: 
are the directions of the easy-axis and easy-plane respec- 
tively, and and are the strengths of the easy-axis 
and easy-plane anisotropy fields respectively. The hight 
of the energy barrier between the two easy-axis directions 
is Eb = tioH,M,s/2. ^ 

The spin-current Is{t) consists of two pulses. An AC 
pulse if^' sin{ujt), followed immediately by a DC pulse 
if-^ as shown in Fig. [2j These pulses are character- 
ized by six parameters: the AC driving frequency w, the 
AC and DC pulse durations, tAc and tjjc, the AC and 
DC spin-current amplitudes, and if *^ (in units of 
magnetization) , and the polarization vector nip such that 

Is — Is nip. A sample switching trajectory is shown in 
Fig. [2] for a free layer with uniaxial anisotropy. 
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FIG. 2: (Color online) Simulated switching trajectory cal- 
culated via numeric integration of the LLG equation; 
(black, solid) and energy E (blue, dotted) along with spin- 
current in arbitrary units (red,dashed). Here time is in units 
of {'yMs)~\ If'^ = 1°^ = 0.02M, and lj = -yH^, with other 
parameters the same as in Fig. [1] 

As a specific example we performed numerical simula- 
tions of a free layer with uniaxial anisotropy, i?^ = 0.5Ms 
and Hx = (taken from Ref. T5] ) , and polarization vec- 
tor nUp — —Gx — (for schematic of MTJ see inset Fig. 
[1]) . For simplicity the AC and DC spin-current strengths 
are taken to be equal Is = Xf~^ = ■ The negative 
sign in the polarization vector is to ensure switching from 
the to —Gz energy minimum using positive DC cur- 
rent. 

To determine the optimal switching protocol numer- 
ous trials were done for a range of spin-current ampli- 
tudes, frequencies, and AC/DC pulse times. Each trial 



consisted of 10^ simulated switching events. The en- 
ergy lost to Joule heat is J = /q*°'°' R{ih)X'^ {t)dt where 
ttotai = tAC+toc is the total pulse time and i?(rh) is a re- 
sistance that depends on the magnetization direction. In 
practice this resistance varies little through the switching 
process [Ml, [13 thus the Joule heat may be approximated 
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The parameter set resulting in the least Joule heat loss, 
with switching probability above 99.5%, for each spin- 
current amplitude I^ is shown in Fig. [TJ Notice the 
optimal energy loss using this AC/DC method is almost 
a third of the purely DC method for the same device and 
the optimal current is roughly half. 

To understand this optimal AC/DC protocol it is first 
convenient to describe the LLG equation, in terms of the 
locally orthogonal coordinates of energy E and angle (f 
For spin-currents that are not too large (i.e. on 
order of the DC critical switching current) i? is a slow 
variable relative to see Fig. [21 This separation of time 
scales allows the equation of motion for E to be averaged 
over each StonerWohlfarth (SW) orbit [ll,[2l,[23. For 
DC spin-currents this reduces the LLG equation in the 
initial potential well to 



E 



~aU{E)+IlVDc[E), 



(3) 



where a is the dimensionless Gilbert damping constant 
and If = if^mp • (— e^) is the portion of the DC spin- 
current polarized along the easy-axis, — direction. For 
the uniaxial case the generalized forces, averaged over the 
SW orbit with energy E^ are given by 



U{E) 
Vdc{E) 



l^E{Eb 
-f2E{l - 



- E)/Ms , 
E/Ebf^ 



(4) 



and correspond to the the effects of damping and DC ST 
respectively. 

For AC spin-currents, with driving frequency uj close to 
the natural frequency of the free layer, the magnetization 
tends to precess with the driving frequency uj [26l . [27| . 
As a result the precessional dynamics can be expressed 
in terms of a relative phase (j){t) = ip{t) — uit which is 
also slow varying relative to ip. Thus the same averaging 
procedure used in the DC case may be applied to both 
E and (j) equations of motion [l^, [3, along with the 
relation sm{ujt) = sm{ip — (f) to single out terms even 
and odd in (p, giving 



E = -aUiE) -If Vac (E) sine 
^ = n{E)-uj-IfW{E)cos(l: 



(5) 



in the initial potential well. Here ^{E) is the natural 
frequency of the free layer about the SW orbit with en- 
ergy E, and I^ = I^'-^mp • (— e^,) is the portion of the 
AC spin-current polarized along the easy-plane axis, —Bx 
direction. 



3 




-Till 7t/2 

<^ 

FIG. 3: (Color online) Simulated energy E vs phase <j) 
(black, dashed) calculated via numerical integration of the 
LLG equation, average energy and phase (red, dot-dashed) 
calculated from integration of Eqn. ((5|, and % = Q trajectory 
(blue, solid) from Eqn. ([7]). Same parameters as Fig. [S] 



A sample trajectory for zero initial energy, 'H(i?, 4>) = 0, 
is shown in Fig. [31 Notice the simulated trajectory from 
the LLG equation and the trajectory calculated from Eqs. 
([S]) closely follow the H — trajectory making it a good 
zeroth order approximation for weak damping, a <C 1. 
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For the uniaxial case fl{E) — ^Hz (1 
the two new generalized forces are given by 



EjEbf^'^ and 



Vac[E) 



j{E-Eb){E/Et)'^^ 
-f{Et/4Ef' , 



(6) 



which correspond to the effects of the AC ST on the 
energy and precessional motion respectively. A sample 
trajectory for E and </> simulated using the LLG equation 
is shown in Fig. [3]along with the trajectory given by Eqs. 
([5]). Notice the energy first overshoots the equilibrium 
energy before winding down to it. 

Setting the LHS of Eqs. ([5]) to zero give the equilib- 
rium energy and phase. The equilibrium energy E^q for 
the uniaxial case is shown in Fig. 2] along with simulated 
values. The hysteresis curve formed is in good qualita- 
tive agreement with experimental and macrospin results 
in Ref. [26l [28l |. The simulations show two jump fre- 
quencies at wi and uj2 where the energy abruptly jumps 
between the upper and lower branches. Since the initial 
energy is near zero for any switching process, for lj < ijJ2 
the trajectory will overshoot and then relax to the lower 
equilibrium branch while for lo > UJ2 the trajectory will 
overshoot and then relax to the upper. As a result a;2 ^ a; 
produces the largest overshoot energies. 

In the absence of damping, a = 0, the trajectories 
given by Eqs. (O possess an integral of motion and can 
thus be described as lines of constant value for some func- 
tion ?^(£:, (/.) For the uniaxial case 
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FIG. 4: (Color online) Equilibrium energy as a function of 
frequency found via numeric integration of the LLG equation 
(red,triangles) and calculated from Eqn. ([Sjl. Same parame- 
ters as in Fig. [2l 
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FIG. 5: (Color onhne) Plots Eac (red,dashed) and 2Edc 
(black, solid) in arbitrary units for parameters = Ms , = 
0.25Ms, a = 0.015, = X| = 0.02Ms, nip = -e^ - and 

Using Eqs. P^ . ^ and (O the energy trajectory for 
any AC/DC spin-current can be calculated. However, the 
question of when to go from an AC to a DC spin-current 
is still of critical importance. To solve this we look for the 
combination of AC and DC pulse times witch minimize 
the Joule heat loss for a given set of parameters. These 
pulse times may be calculated using Eqs. ([T6| and ^ as 
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Eac 
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Be EdC 



(8) 
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where the (p dependence in Eac is eUminated using Eqn. 
d?]), and Ec is the energy where the spin-current goes 
from AC to DC. In this description the optimal pulse 



times are found by solving = 0. Substituting Eqs. 
([8]) into Eqn. ([2]) and taking the derivative with respect 
to Ef. gives the optimization criteria 



Edc{Ec) _ ^Eac{Ec 



(9) 



Thus the optimal strategy when if^ — I^^ is to use 
an AC spin-current until the DC energy current is twice 
the AC energy current. An example of this realization is 
shown in Fig. [5j It should be noted that Eqn. © may 
have more than one solution. Thus care must be taken 
to ensure the E,. used is indeed the global minimum. 



*AC^*total 




AC/DC method Fig. |6] shows the optimal AC pulse 
time as a fraction of the total pulse time tAc /Uotai , for 
a broad range of anisotropy configurations H^/Hz and 
spin-current polarizations ratios I^/If . Here tAc/Uotai 
is calculated numerically using Eqs. ([9]) ([8]) for AC fre- 
quency oj — n{0) and Xf-^ — I^^ — 2Xc where Ic is the 
DC critical current. One can see there is a large region 
where purely AC spin-current pulses are more energy ef- 
ficient than purely DC pulses {tAclUotai = 1), as well as 
a broad a broad region where the optimal strategy uses 
combined AC/DC spin-currents, (0 < tAc/Uotai < !)■ 
There is also a narrow range where DC spin-currents are 
optimal {tAc/ttotai — 0). This range of results stresses 
there is no single general optimal application of either 
AC, DC, or combined AC/DC current pulses for MTJs. 
Instead the optimal strategy for each MTJ must be de- 
termined carefully as outlined above. 

In conclusion we have developed a theoretical descrip- 
tion for magnetic switching using consecutive AC and 
DC spin-current pulses as well as a means for minimiz- 
ing the Joule heat loss associated with such strategies. 
It provides a means of determining the optimal AC and 
DC pulse durations, tAc and t^c, as well as the op- 
timal AC driving frequency uj2 ^ ^opt and spin-current 
strength. Such optimized strategies may use currents well 
below the DC critical current and result in significantly 
less Joule heat loss than purely DC or AC spin-current 
methods. 

This work was supported by NSF Grant DMR-0804266 
and U.S. -Israel Binational Science Foundation Grant 
2008075. 
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FIG. 6: (Color online) Shows tAc/ttotai for a range of H^/ 
and FJII. Here H^, = 0.5M,, a = 0.015, u = 0(0) and 
= 2Xc where is the DC critical current. 

Using Eqs. ([8|) ([9|) the Joule heat loss can be cal- 
culated as a function of spin-current and frequency 
J(I^*-^,lf^*-^, oj). The frequency w that gives most effi- 
cient switching protocol does not lend itself to a simple 
analytic solution, however determining it numerically is 
relatively simple. The Joule heat loss using this optimal 
frequency is shown in Fig. [T]as a function of spin-current 
along with simulations of the LLC equation, optimized 
numerically for a uniaxial free layer. From simulations 
and the analysis described here, we found the optimal 
frequency to be slightly larger than the upper bifurca- 
tion frequency uj2. 

The optimization method for this AC/DC method out- 
lined above is not limited to the uniaxial case shown 
here. It can be applied much more generally allowing 
an optimal AC/DC spin-current protocol to be deter- 
mined for any free layer with easy-axis and easy-plane 
anisotropy, see supplemental material [l6j . To illustrate 
the range of possible strategies that may optimize this 



The motion of the magnetic unit vector rh{t) is de- 
scribed by the Landau-Lifshitz Gilbert equation with 
Slonczewski spin torque term 



dm 
~dt 
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(10) 



where 



Tgd -7am X [m X Hcff] , 
Tsion = -7X5 (i)m X [m X nip] 



(11) 



are the conservative Landau-Lifshitz torque, Gilbert 
damping torque, and Slonsnewski torque respectively. 
The magnetization can also be described by the coor- 
dinates energy E and angle Lp. The angle is the time 
into the orbit with energy E in the absence of spin-torque 
and damping, normalized to 27r and calculated from 



dt = 



■ LL 



dm 



■ LL 



dip 

n{E) 



(12) 



where Q{E) is is the precessional frequency of the orbit 
with energy E such that § dt = 2-11 /Q,{E). Throughout 
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this appendix Lp is taken to be zero along the easy-plane 
axis, e^,. 

The magnetization rh{E, ip) is related to E and (p by 
the partial derivatives 



dEiaci = 7— ^ . (13) 



■ LL| 



The unit vectors E and tp are locally orthogonal, d,^m ■ 
dEi&i- = 0, thus the LLG equation may be written as 
dxh/dt ~ d^ihip + Betx^E. Separating our the energy 
and angle equations of motion gives 



E — — Hoff ■ (Fgd + Tsion) ; 

TlL • Tsion 



^ ^ n{E) - n{E)- 



(14) 
(15) 



LL 



Averaging Eqs. p4)) and p5|) about p for DC spin- 
currents gives 



E = -aU{E)+I,VDc{E), 

ip = n{E) 



(16) 



where 



UiE) 



n{E) 



j> [dm X Heff] • M , 

Vdc{E) = j> [dm X M] • mp . 



(17) 



Here the contour integrals follow orbits of constant en- 
ergy. 

For AC spin currents this averaging gives 



E = -aU{E) -IsVAc{E)smc 
(j) = n{E) -uj -I,W{E) cos 4 



(18) 



where (j) = ip — uj and 

Vac{E) ^ [dm cos (19) 

mE)=l^§ (-'^■^l-^^H^^^-.) sin p . (20) 

Here again the contour integrals follow orbits of constant 
energy. 

In the absence of damping a = Eqs. p8|) follow lines 
of constant value (j)) given by 



■H[E, 0) = y dE'j{E') \uj - n{E') + 1, W{E') cos < 


(21) 

where J{E) is the solution to the linear differential equa- 
tion 
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